PROCEEDINGS 

™ SCIENCE 



Exploring the epsilon regime witli lattice Wilson 
fermions 



o 
o 

(N 
u 

< 
00 



i3 



> 
O 

o 

O 

> 
• ^ 

X 
S3 



Oliver Bar", Silvia Necco*' , Stefan Schaefer" 

" Institute of Physics, Humboldt University Berlin 

Newtonstrasse 15, 12489 Berlin, Germany 

^ CERN, Physics Departement, 1211 Geneva 23, Switzerland 

E-mail: obaer@physik.hu-berlin.de, Silvia . NeccoScern . ch, 
sschaef Sphysik . hu-berlin . de 

We study the impact of explicit chiral symmetry breaking of lattice Wilson fermions on mesonic 
correlators in the e-regime using Wilson chiral perturbation theory. We generalize the e- 
expansion of continuum chiral perturbation theory to nonzero lattice spacing a and distinguish 
various regimes. It turnes out that lattice corrections are highly suppressed, as long as quark 
masses are of the order aAQQ^. The lattice spacing effects become more pronounced for smaller 
quark masses and may lead to non-trivial corrections of the continuum results at next-to-leading 
order. We compute these corrections for standard current and density correlation functions. A fit 
to lattice data shows that these corrections are small, as expected. 
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1. Introduction 

Lattice QCD simulations with light quarks are now approaching domains where a reliable 
matching with the chiral effective theory can be performed. Through the matching it is possible 
to extract the Low-Energy Couplings (LLCs) of the effective theory: many results have been pre- 
sented in the past months for the leading order (LO) and next-to-leading order (NLO) couplings, 
both for the Nf = 2 and Nf = 3 effective theory (for a recent review see for instance [1]). A very 
important issue which concerns these determinations is the control over the systematic uncertain- 
ties. From this point of view, it is very useful to extract the LECs from a large set of observables 
and from different kinematic regimes: this will give solidity to lattice results and will help to get a 
comprehensive picture of low energy properties of QCD. 

An interesting approach is to investigate QCD in a finite volume V = L^T in the so-called e- 
regime [2, 3], where the pion wavelength is larger than the size of the box, Mj^L < 1. The relevant 
feature of this regime is that, due to a reorganization of the chiral expansion, volume effects are 
enhanced, while mass-effects are suppressed with respect to the usual infinite-volume case (or p- 
regime, where Mj^L » 1). For this reason, at a given order in the perturbative expansion, less LECs 
will appear: predictions are less "contaminated" by higher order unknown couplings, making the 
£-regime potentially convenient for the extraction of the LO constants. 

From the point of view of lattice computations, reaching the £-regime requires simulating 
small quark masses; this may influence the choice of the discretized Dirac operator to be adopted. 
Ideally, a Dirac operator which satisfies the Ginsparg- Wilson relation [4] would be very advanta- 
geous, since it preserves chiral symmetry at finite lattice spacing, and small quark masses are acces- 
sible. Moreover, the topological charge can be unambiguously defined through the index theorem, 
and the lattice results may be matched with the predictions of the chiral effective theory at fixed 
topology [5]. A large number of quenched computations in the £-regime with Ginsparg -Wilson 
Dirac operator has been performed (see [6] for a recent study, and references therein for precedent 
computations). Even though the results were to a large extend promising, the main obstacle for 
progress in real QCD is the fact that simulations with dynamical sea quai^ks ai^e still extremely 
time-consuming. For recent calculations with dynamical chiral fermions see [7, 8, 9, 10]. 

On the other hand, dynamical simulations with the Wilson Dirac operator with 0{a) improve- 
ment are becoming fairly inexpensive. Reaching small quark masses with Wilson fermions has 
been considered problematic for many years; this issue can be faced by adopting techniques such 
as reweighting [11]. This method has been applied in [12] allowing to reach the £-regime. Also the 
ETM collaboration investigated the £-regime with a twisted mass Wilson Dirac operator [13]. 

Since Wilson fermions explicitly break the chiral symmetry at finite lattice spacing, the match- 
ing with the chiral effective theory should be performed only after a continuum extrapolation of 
the lattice results. While this is not an unrealistic goal for the near future, the presently available 
simulations in the £-regime are canied out at a single value of the lattice spacing. In [12] the pseu- 
doscalar and axial coiTclations functions turned out to be very well described by the predictions of 
the continuum chiral effective theory at NLO. Similar observations have been made by ETM. Still, 
it is important to have a theoretical understanding of the impact of explicit breaking of chiral sym- 
metry on computations in the £-regime. We address this question in [14]: the tool that we adopt 
is the so-called Wilson Chiral Perturbation Theory (WChPT) [15, 16], the low-energy effective 
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theory for lattice QCD with Wilson Dirac operator. A similar analysis has been carried out in [17]. 



2. Wilson Chiral Perturbation Theory 

The chiral effective Lagrangian of WChPT is expanded in powers of pion momenta p^, the 
quark mass m and lattice spacing a. Based on symmetries of the underlying Symanzik action 
[18], the chiral Lagrangian including all terms of 0{p^ , p^m,rrP' , p^a,ma) is given in [16]. The 
0{a^) contributions are constructed in [19, 20]. In the following we will restrict ourselves to the 
case Nf = 2 with degenerate quark mass ni. The leading order Euclidean chiral Lagrangian in the 
continuum is given by [21, 22] 

if2 = ^Tr (5^ [/^^ [/^) - ^— Tr + 1^). (2.1) 

The pseudo Nambu-Goldstone modes are parametrized as usual by the SU(2) field 

[/(^)=exp(2/^(x))/F, 

and F , B are the familiar LO couplings. The leading terms involving the lattice spacing are 

= aW45Tr {df^UdfiU^) Tr {U + U'^) - dmW^s (Tr {U + U"^))^, (2.2) 
^,2 = ^C2«2(Tr([/ + [/t))', (2.3) 

where m = 2Bm and d = IWqo. W45, Wgs, Wq and C2 are new LECs which are not determined by 
the symmetries. Note that the mass parameter m in Eq. (2.3) is the so-called shifted mass [15]: 
besides the dominant additive mass renormalization proportional to 1 /t? it also contains the leading 
connection of 0{a). 

Currents and densities in WChPT can be constructed by a standard spurion analysis or by 
introducing source terms. Here we report the axial vector cuiTcnt and the pseudoscalar density 
including the leading 0{a) corrections [23, 24]: 

^M,wChPT = ^^,cont { 1 + [W45Tr([/ + [/^) +4Wa] | + 2dWiod^,Tr {T%U - [/+)) , (2.4) 

^ChPT = P"oni { 1 + [W(>{Tv{U + U'^) + AWp] I , (2.5) 

where 

^M,cont = yTr [T'iU^d^U - Ud,U^)) , P« „t = '^Tr (r ([/ - f/t)) , (2.6) 

and T" are SU(2) generators normalized such that Tr(rT^) = 5''*/2. Notice that the LECs Wa,p 
stem from the renormalization factors, which up to 0{a) take the form p = 1 + l6dWA.p/F^. 
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2.1 Power counting in infinite volume 

In WChPT there are two parameters which break explicitly the chiral symmetry, the quark 
mass m (counted as 0{p^)) and the lattice spacing a. The power counting is determined by the 
relative size of these two parameters. In particular, one distinguishes [23, 25] two different regimes: 
(i) the GSM ^ regime, where m ~ fl'A^f.Q and (ii) the Aoki regime where m ~ o^A^q^^. In the Aoki 
regime lattice artefacts are more pronounced, and the in Eq. (2.3) enters already at LO. The 
pion mass at leading order is given by 

GSM regime : Ml = 2Bm, (2.7) 
Aoki regime : Mq = 2Bm - 2c2a^ ■ (2.8) 

The sign of C2 governs the phase diagram of the theory. The reader can refer to [15] for a complete 
discussion. 

2.2 Power counting in the e -regime 

As already anticipated in the introduction, the £-regime is a finite-volume regime where the 
pion wavelength is larger than the size of the box, M^L < 1 (but still L » 1 /Aqcd) [2, 3]. This 
con^esponds to approaching the chiral limit by keeping the dimensionless quantity pL = mXV < 
0(1) (where £ = F^B is the quark condensate in the chiral limit). The main effect of formulating 
the effective theory in this regime is that the pion zero-mode becomes non-perturbative and its 
contribution has to be treated exactly. This is achieved by factorizing the pseudo Nambu-Goldstone 
boson fields as 

[/(x) = exp(^|^(x)^[/o, (2.9) 

where the constant Uq € SU(2) represents the collective zero-mode. The non-zero modes t, can 
be still treated perturbatively. The £-regime requires a reorganization of the perturbative series: in 
the continuum, this corresponds to taking the quark mass of order m ~ 0(£^). Mass effects are 
hence suppressed compared to the p-regime (or infinite volume) case, while finite-volume effects 
are enhanced and become polynomial in (FL)^^. 

We now want to extend the WChPT to the £-regime. Also in this case we have to assign a 
relative power counting of the lattice spacing a with respect to the quark mass m. If we assume that 
the quark mass can be considered of order m ~ 0(£4) also in WChPT 2, we obtain 

GSM regime : m ~ 0[aA\cv>) ar^ 0(£'^), (2.10) 
Aoki regime : m 0{a^ I^qq-q) a^O{e^). (2.11) 

Moreover, the £- expansion allows us to introduce another intermediate regime between the GSM 
and the Aoki regime: we can define the GSM* regime, where a ~ 0{£^). 

We are interested in computing two-point correlation functions within the WChPT in the £- 
regime. In particular, we give explicit results for the pseudoscalar and axial time correlators, 

5"'Cpp{t) = jd'x{P"{x)P\Q)), 8"'CAA{t) = |^3x(Ag(x)Ag(0)). (2.12) 

'GSM stands for generically small masses. 

^While this is a natural choice in the GSM regime, the situation in the Aoki regime can be more subtle. See [14] for 
a more detailed discussion on this subject. 
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Currents and densities are defined in Eqs. (2.4, 2.5); the subscript "WChPT" is now omitted. By 
adopting the power counting that we have defined, we find that in the GSM regime lattice coitcc- 
tions enter only at NNLO, while in the GSM* regime they appear already at NLO. In particular, the 
leading connection is given only by the term ^^,2 in Eq. (2.3). The reason of this suppression can be 
traced back to the fact that the lattice spacing conections in the chiral effective theory action and 
in the effective operators are either quadratic in a or they come with an additional power of either 
m or p^. Hence, in the e-regime the suppression of lattice spacing corrections works similarly as 
the suppression of mass effects. On the other hand, in the Aoki regime effects of lattice artefacts 
are more severe and show up already at LO. 

It is important to observe that these considerations are valid for unimproved Wilson fermions. 
If the theory is non-perturbatively 0{a) improved, the corrections due to as well as the 0{a) 
terms in the operators are absent, and lattice artefacts are due to only. Consequently, in the e- 
regime the leading corrections due to lattice artefacts are essentially unaltered for the unimproved 
theory, since improvement acts only on subleading terms. 

2.3 Leading corrections in the GSM* regime 

The continuum pseudoscalar and axial correlators at NLO in the e-expansion can be written 
as [26] 

Cpp^AA ct(0 = ap^A + bp^Ahi (t/T) , (2. 13) 



where 
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"1-2 



1 

12 



(2.14) 



For Nf = 2 the coefficients ap^A,bpA explicitly read [26] 
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I\, I2 are modified Bessel functions of the first kind; /3i and ^00 ai"e so-called shape factors [27, 26], 
which depend only on the geometry of the finite box. Zgff is the quark condensate at one loop [3] 



Teff = r 1 + 



2F^JvJ' 



(2.17) 



and Hen = niLsffV. As already anticipated, the continuum NLO predictions in the e-regime contain 
only the LO LLCs £ and F. 

The first non-trivial modification of the continuum NLO results appear in the GSM* regime 
and it is due to only. In this case we can write down the full NLO correlators in WChPT as 



Cpp^AA{t) — Cpp^AAct{t) +CppAAcfi{t). 



(2.18) 
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By performing the explicit computation (see [14] for the full details) it turns out that the corrections 
CppAAa-{f) ^1'^ time-independent and hence affect only the constant part of the con^elators. In 
particular we obtain: 

CpPcfl{t) = —p\al, CAAa^{t) = ypA,2, (2.19) 

where 

^ ^ 5^7? (2m) - 10/i(2m)/2(2m) - 3^/1(2^) 
2At3/2(2At) 

and p = F^C2a^V is the dimensionless LEC which parametrizes the 0{a^) correction. 
It is useful to compute the leading 0{a^) corrections to the PCAC quark mass: 



mpcAC = m 



1+P 



2 /i(2m) \ 
At2 lihilli)) 



(2.21) 



It is now possible to express the correlators Cpp^AA{t) as a function of /i = mpcAC^^; the result is 

Cpp{t) = Cppctit) + ^pA„2, Caa(0 = CAActit) + ypAa^-, (2.22) 

where the continuum correlators are as in Eq. (2.13), but with the replacements ji ^ fx, /leff — > 
Aeff = nipcAC^eff^ and 

Other correlation functions can be computed along the same line. For instance, in [14] we report 
also the result for the vector correlator. 

3. Reanalysis of lattice data and conclusions 

These predictions from WChPT at NLO can be tested against lattice data generated in [12], 
where pseudoscalar and axial correlators have been computed on an ensemble with Nf = 2 flavours 
of dynamical improved NHYP Wilson fermions [28]. The lattice spacing is a ~ 0.1 15 fm, and two 
lattice extents are available, Li = 16(3 ~ 1 .84 fm and Lj = 24^ ~ 2.8 fm. Quark masses approach the 
£-regime, with jx ~ 0.7 — 2.9 for the volume Vi = L\ and /i ~ 2.1 — 5.0 for the volume V2 = L\. In 
the GSM* regime, we have only the additional LEC C2 with respect to the continuum case. Notice 
that its value will depend on the particular discretized action which is used. We simultaneously fit 
the two correlators for all available quark masses; for the volume V2, a fit in the range t G [6, 18] 
gives 

i:^(At = 2 GeV)j = 249(4) MeV, = 88(3) MeV, C2 = 0.02(8) GeV^. (3.1) 

The data, along with the theoretical curves, are shown in Fig. 1. The errors from the renormal- 
ization factors Za, Z^^(/x = 2 GeV) computed in [12] are not included in the uncertainties of the 
LECs. Varying the time range of the fit does not give significant differences for the LECs, as long 
as fmin/fl' > 4. Also discarding the heaviest mass does not change the results of Eq. (3.1) within 
the statistical errors. The smallest volume V\ yields values which are consistent with Eq. (3.1), but 
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Figure 1: Fit of the WChPT predictions to lattice data. All data points within the 
fit range of t/a E [6,18] for the four sea quark masses are included in the combined 
fit. The hopping parameter K = (0.128150,0.128125,0.1281,0.128050) corresponds to ampcAC = 
(0.0019(4), 0.0024(3), 0.0030(3), 0.0044(3)) respectively. The axial vector coiTelator is multiplied by a 
factor 50 for better visibility. 

the large of the fit may indicate that NLO formulae are no longer applicable. The values of F 
and £ are compatible with other determinations [1], while the value of C2 is compatible with zero. 
A continuum fit (with C2 = 0) yields virtually unchanged values for F and £, showing that cut-off 
effects do not impact the extraction of the LECs beyond the level of the statistical uncertainties. 

This is a very encouraging result: simulations with Wilson fermions in the £-regime are fea- 
sible and seem to be a viable alternative to dynamical simulations with chiral fermions. Similar- 
conclusions have been reached in [17]. The results derived here can be generalized in various ways, 
for example to the case with a twisted mass term or to an ai^bitrary number of flavors. 
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